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In a classical, quartic field theory with SU (N) x Z2 symmetry, a class of kink solutions can be found 
analytically for one special choice of parameters. We construct these solutions and determine their 
energies. In the limit N — > 00, the energy of the kink is equal to that of a kink in a Z2 model 
with the same mass parameter and quartic coupling (coefficient of Tr($ 4 )). We prove the stability 
of the solutions to small perturbations but global stability remains unproven. We then argue that 
the continuum of choices for the boundary conditions leads to a whole space of kink solutions. The 
kinks in this space occur in classes that are determined by the chosen boundary conditions. Each 
class is described by the coset space H/I where H is the unbroken symmetry group and / is the 
symmetry group that leaves the kink solution invariant. 



I. INTRODUCTION 

Classical solutions in field theories have been of long- 
standing interest and have played an important role in 
our understanding of a variety of phenomena. The exis- 
tence of solutions can be predicted on the basis of topol- 
ogy though the actual construction of the solutions is 
usually quite difficult and requires a certain amount of 
guesswork. Even after a solution is found, one needs to 
check its perturbative stability. And even after a solu- 
tion is proved to be perturbatively stable, there is no 
guarantee that it is the globally stable solution. 

In this paper, we shall construct a class of kink so- 
lutions in purely scalar, SU(N) x Z2 field theories. The 
analytic construction is made possible by a special choice 
of parameters. We then prove the perturbative stability 
of the kink solutions. As there is no known extension 
of Bogomolnyi's method fj|] to this case, global stability 
will remain unproven. We will then address the question 
of whether there might be other kink solutions in the 
model. We discuss the possibility that various bound- 
ary conditions imposed at spatial infinity, while topolog- 
ically equivalent, lead to different kink solutions. Each 
kink solution is continuously degenerate with a degener- 
acy determined by the invariance group of the solution. 
We describe the problem of mapping the space of kink 
solutions. 

The following few sections deal with the construction 
and perturbative stability of the kink solutions that have 
been found. The space of kinks is described in Sec. VII 



and this can be read with minimal reference to the prior 
sections. 



II. THE MODEL 



We will consider the 1+1-dimensional, classical field 
theoryf] 



L = TrP^) 2 ] - V{$) 



(1) 



V(<f>) = -m 2 Tr[$ 2 ] + h(Tr[<^ 2 }) 2 + ATr[$ 4 ] + V (2) 

where $ transforms in the adjoint representation of 
SU(N) (N = 2n + 1) and hence is an N x N hermitian 
matrix. <I> can also be written in terms of components: 



jv -1 

£ 

a=l 



(3) 



where T a are the (traceless) generators of SU(N) and 
will be taken to be in the Gell-Mann representation j|] 
and normalized so that 



TrpyTft] = -S ab 



(4) 



The field components $ a are real. The constant Vo is 
chosen so that the minimum of the potential is at V = 0. 
Note the absence of a cubic term in the model. 

The center of SU(N) is Zjv and the group elements 
of the center transform $ by multiplication by factors 
of exp(i2ir/N). For even N, the transformation $ — > 
— $ is included in the center of SU(N). However, for 



*The addition of gauge fields can easily be accommodated 
in what follows as discussed in Sec. VI 
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odd N it is not, and the full symmetry of the model is 
SU(N) x Z 2 . We will only be interested in odd values of 
N since it is this case that the model admits topological 
kink solutions. 

We will consider parameters h and A such that the 
symmetry breaking is: 

<™ nv7 SU(n+l)xSU(n)xU(l) 
SU(N)xZ 2 ^ - — - (5) 

where N = 2n + 1 (n is a positive integer). This symme- 
try breaking pattern is achieved if [|| 



> 



N 2 + 3 



X N(N 2 -l) 
In fact, in the next section we will choose 
h 3 



A N(N-l) 



(6) 



(7) 



This choice is consistent with eq. (g) as long as N > 3. 
To fix the constant Vo, let the minimum of V(<I>) occur 

at 



$ = $ = R diag(pi,p 2 , -,Pn) 



(8) 



where R is an overall normalization factor and the matrix 
is normalized so that 



N 

E 



i 



pi 



Then the potential is extremized for 
R 

where 



m 

7% 



N 



Then, V($o) = gives 



N 2 



N(N 2 - 1) 



-A . 



_ m 



(9) 



(10) 



(11) 



(12) 



The condition that A' > is precisely the condition in 
eq. (|) 

We shall assume the parameters as constrained by eq. 
( 6 ) , in which case the unbroken symmetry is given in eq. 
(5) and the vacuum expectation value of $ is 



$o = R\ 



nl n+ i 
N(N 2 - 1) V 





-(n + l)l„ 



(13) 



where 1„ denotes the n x n unit matrix, N = 2n + 1, i? 
is defined in eq. ([l0|). 



III. KINK SOLUTION 

The model contains the breaking of a discrete Z 2 sym- 
metry whenever N is odd. Hence kink solutions ex- 
ist. Across these kink solutions, the vacuum expectation 
value of $ must be related by an element of the Z 2 group. 
Hence, if $^ denotes the kink solution, 



- DC 



) = -U<f> k (x = +(X))U- 1 



where U is an element of SU (N) . 

Based on the experience with kinks in SU(5) we 
conjecture the following form for the kink solution: 



$ k (x) = f{x)M + g(x)P 



(14) 



where f(x) and g(x) are unspecified functions as yet, and 
M and P are SU (N) generators such that 



Tr(M 2 ) = - = Tr(P 2 ) 



Tr(MP) = = Tr(M 3 P) = Tr(MP 3 ) . 
Insertion of $^ in the potential gives 



(15) 



(16) 



A[/ 4 Tr(M 4 ) + .g 4 Tr(P 4 )] + 
[^+6ATr(M 2 P 2 )]/V-^o 



(17) 



The important realization that permits an analytic solu- 
tion to be found is that the cross-term containing both 
/ and g disappears if we choose 



-12ATr(M 2 P 2 ^ 



(18) 



In this case, the energy of the kink will separate into two 
pieces, one depending only on / and the other depending 
only on g. So the energy will be a sum of two single 
field energies, each of which can be treated directly or by 
Bogomolnyi's method. 

Let us further choose the matrices M and P in such a 
way that the kink boundary conditions imply 

/(+oo) - -/(-oo) , 9 (+oo) = +.g(-oo) . (19) 

Then, we find that / and g must satisfy 

/" + m 2 / + 4A[3Tr(M 2 P 2 ) - Tr(M 4 )]/ 3 = (20) 



1 



ff 



4A 3Tr(M 2 P 2 ) - Tr(P 4 ) 
The equation for / will have the solution 



f{x) = /otanh (j^= x 



(21) 



(22) 



2 



where 



/o 2 = 



1 



4A Tr(M 4 ) - 3Tr(M 2 P 2 ) ' 
The solution for / is valid provided 

Tr(M 4 ) > 3Tr(M 2 P 2 ) 
and the (constant) solution for g is valid provided 
Tr(P 4 ) > 3Tr(M 2 P 2 ) 



(23) 



(24) 



(25) 



In addition to the properties in eqs. ( |l5| ) and ([l6]), the 
conditions in eqs. (|24|) and (25), M and P must yield 
the desired vacuum expectation values for $ at spatial 
infinity. A choice of M and P that satisfies all these 
conditions is: 



M= (3 




where, 



y/2N(N-l) 



(26) 



(27) 



Straightforward calculation shows that all the conditions 
on M and P are met provided that N > 3. 

With this choice of M and P, the kink solution can 
now be written as: 



m 

71 



N -1 



N 



tanh I —=x M 
V2 



NP 



(28) 



An alternative, more transparent, form of this solution 
is: 



$fc(x) 
where 



f - F(x) 



1 + F(x) 



1', 



F{x) 



tanh | — =x 

,\/2 . 



(29) 



(30) 



and <I>± = $(x = ±oo). Note that the alternate form 
does not work for any chosen boundary condition. For 
example, if $+ = — <£>_, the form leads to the embedded 
kink solution which is known to be unstable Q] . 

Now that we have shown that <£>fc is a solution within 
the restricted ansatz in eq. (|l4|), we also need to show 
that it is a solution of the full theory. This is most simply 
done by writing 



$ = $ fc + * = $ fe 



(31) 



and then checking that the energy density does not con- 
tain any terms that are linear in ^ . In the quadratic 
terms in the energy density, this will clearly by the case 
since satisfies the equations of motion ( p0| ) and ( pl| ) 
and the generators satisfy the orthogonality condition in 
eq. ([|). The only terms that can potentially lead to a 
term linear in the ip a come from the quartic term, Tr($ 4 ), 
in the potential and are of the form 

Tr(<^T Q )V> a 

For off-diagonal T a , this vanishes since $| is diagonal 
and a product of a diagonal and an off-diagonal matrix 
is off-diagonal. For diagonal T a , it is better to choose 
a representation of the generators other than the Gell- 
Mann representation. The choice of diagonal generators 
other than M and P are: 



V 



T 













(: 


o„ 





























-1 





) 













(t = l, 



1) (32) 



,2n-2) 



(33) 



where Tj (i — l,..,n— I) are the normalized, diagonal 
SU{n) generators in the Gell-Mann representation. The 
off-diagonal generators are the same as in the Gell-Mann 
basis. It is easy to check that this basis is complete and 
the generators satisfy the normalization in eq. (Q) . 

Now we are interested in checking if there are terms 
in the potential that are linear in the components of ^ 
expanded in the new basis. It is easy to check that 



Tr(M 3 T ! ) = = Tr(P 3 :P) 



Hence 



Tr($|T !; 







and there are no linear terms in the components of ^ 
occurring in the energy density. The terms linear in the 
perturbations along the generators M and P will vanish 
because these have already been chosen to satisfy the 
equations of motion. 

Therefore there are no linear terms in ^ in the energy 
density and in eq. (p8|) is indeed a solution. 



IV. KINK ENERGY 



The energy of the kink is: 

E k = J cLT[Tr(<i>' fc ) 2 + V($ k )} (34) 
Insertion of $fc from eq. (Eq) and evaluation yields 



3 



2^2 m 3 



N - 1 



3 A \N-3 
In the limit that N — > oo, this gives 



2 V / 2to 3 

~3 A~ 



(35) 



(36) 



which is precisely the energy of the kink in the Z 2 model: 



L 



A 



m 
4A 



2 r ' 4 

In the large iV limit, the kink solution (E8fi goes to 



(37) 




tanh(ma;/\/2) 



(38) 



Note that this is not traceless because we have discarded 
a large number (N) of small (order 1/N) terms. 

V. PERTURB ATIVE STABILITY 

The procedure for proving perturbative stability is 
straightforward though tedious. We consider small de- 
viations \& from the kink solution as in eq. ( |3~l| ) and then 
find the change in the energy due to the perturbations, 



8E[V] = Tr J dxV 



(39) 



where the matrix V2 is obtained by expanding the po- 
tential V up to quadratic order in \I>. If the Schrodinger 
equation 



dx 2 



, , V 2 ($ fe ) 



(40) 



does not have any negative eigenvalues - i.e. regular 
solutions where \P — > at spatial infinity only exist for 
u> > - then the solution is perturbatively stable. 

The tedious part of this calculation is the evaluation of 
V 2 ($fc). Note that V 2 is an (N 2 - 1) x (N 2 - 1) matrix 
since there are iV 2 — 1 components of 'J, one for each 
generator of SU(N). Let us write 



N'-l 

0=1 



(41) 



and the elements of V2 as Viab- 

Let us discuss off-diagonal perturbations first. These 
fall into 5 types. The first type is when 



Pn 

Tl cx I 




(42) 



where p n denotes a non-trivial, off-diagonal nx n matrix 
and the trivial n x n matrix. The second to fifth type 
of perturbations are: 



T? oc I p\ 



T 3 a cx 




(43) 



(44) 



(45) 



(46) 



where c is a non-vanishing complex n-dimensional mul- 
tiplet. Then the elements of V 2 for Type 1 and Type 4 
perturbations are: 



,N±3F 



2aa 



N -3 

and for type 3 (minus signs) and 5 (plus signs) are: 

o 

77? 

V 2aa =±—F(1±F) 



(47) 



(48) 



where the function F(x) is defined in eq. (J30|) - The 
elements of V 2 vanish for type 2 perturbations and there 
is no mixing between different off-diagonal perturbations 
either. 

The potential V<2, aa for type 1 and 4 perturbations 
is everywhere non-negative for N > 3 and hence the 
Schrodinger equation will not have any bound states. 
The situation for type 3 and 5 perturbations is less clear 
since ±F(1 ± F) can have either sign. However the 
Schrodinger equation has the u> = solution tp = 1 =p F 
and this solution has no nodes i.e. "0^0 except at in- 
finity[]. The zero mode solution does not go to zero at 
either x — — oo or at x — +oo. Then, in order for a 
solution to go to zero at both spatial infinities, u> has to 
be chosen to be greater than zero. The reason is that 
the to = solution needs to have smaller curvature (i.e. 
the curvature needs to be more negative) so that it can 
vanish at infinity. The curvature is proportional to —lu 
in the region where the potential is small. Therefore uj 
has to be larger. Hence, once again, there are no bound 



^This zero mode will be of special interest in Sec. VII 
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states to the Schrodinger equation. This shows that the 
solution is stable to off-diagonal perturbations. 

Next we consider diagonal perturbations. These can 
be classified in three types. If we write: 














s„ 














(49) 



where R Q and S a are n x n, diagonal matrices, then the 
three types of T a are: (1) R a ^ 0, S Q = 0, Q a = 0, (2) 
R a ^ 0, S a ^ 0, Q a = 0, and, (3) R Q ^ 0, S a ^ 0, 
Q a 0. Then there are six cases to be treated in finding 
the elements of V2 . 

If both T a and T are type (1), we find 



V 2ab 



, N + 3F 
7V-3 



Sab > 



(50) 



where we note that we are considering TV > 3. Since 
V2ab is everywhere non-negative, the Schrodinger equa- 
tion has no bound states and the solution is stable to 
these perturbations. 

If both T a and T are type (3) - there is only one 
generator of type (3) - we find 



V 2aa = ^-(3F 2 - 1) 



(51) 



Now the Schrodinger equation is precisely that obtained 
when considering fluctuations about a Z 2 kink 0. The 
complete eigenspectrum of this equation is known B and 
the lowest eigenvalue is u> = corresponding to the zero 
mode which describes translations of the kink. So there 
are no bound states and no instability to these perturba- 
tions. 

The only other non-trivial perturbations are when T a 
and T are both type (2). For these perturbations, we 
have 



Ra ^a^-n 



(52) 



where a a is determined from the normalization of T a . 
After some algebra, we find 



^aV 2a b^b 



N-3 



(N-3F)ip 2 ± + {(N-3F) 
60- 2 (N - 1)(1 - F)}Vf ] 



(53) 



where, we have defined the unit vector df,, decomposed 
ijj a parallel (?/>|| = a a ^a) and perpendicular (ip±) to the 
vector a a , and written 



ct 2 = a a a a . 



The coefficient of ip± is positive since N > 3 and hence 
this perturbation cannot cause an instability. The coef- 
ficient of V> 2 needs to be checked. 



Using the normalization of T a , we obtain 

1 



2(iV-l) 



(54) 



Inserting this into eq. (|53|), we find that the coefficient 
of "0 2 is simply +m 2 . Hence there is no instability to 
diagonal type (2) perturbations. 

This explicit analysis shows that the solution is sta- 
ble to all perturbations but this does not imply that the 
solution minimizes the energy globally. 



VI. GAUGE FIELDS 

The inclusion of gauge fields, = A^T a , has no effect 
on the existence of the static solution in eq. ( p8| ) . This 
can be seen by noting that the terms in the energy that 
involve the gauge fields are: Tr(E 2 + B 2 ), -Tr([A x , $] 2 ) 
and iTr(A x [Q, d x &]. The first two terms are quadratic 
in the gauge fields while the last one vanishes because 
d x Q] — for the solution. Hence A^ = is a solution 
of the equations of motion. 

The presence of gauge fields does not have any effect 
on the perturbative stability of the solution. To see this 
one can check that both the quadratic order terms in 
the gauge field are non-negative. For the first term this 
is explicit while for the second term one uses the fact 
[A x , $] = — [A x , $P since A x and $ are Hermitian. Hence 
the kink solution with A^ = is stable to perturbations 
in the gauge fields. 



VII. SPACE OF KINKS 

Different boundary conditions will, in general, lead to 
different kink solutions. Therefore kinks with different 
boundary conditions fall into different classes - kinks be- 
longing to different classes cannot be transformed into 
one another by global SU(N) rotations. Here we would 
like to find the degeneracy of a kink solution i.e. the 
space of boundary conditions that lead to degenerate 
kink solutions. There is clearly an SU(N) global de- 
generacy but this is not very interesting since it applies 
to any field configuration in the theory. It is of greater 
interest to only consider those global SU(N) transfor- 
mations that leave = $(2; = —00) unchanged but 
act non-trivially on $ + = $(x = +00). If we denote 
the unbroken symmetry groups at x = ±00 by H±, such 
transformations belong to But the transformations 
that belong to K = H + n H- will act trivially on both 
$_ and on $ + . Therefore the space of boundary condi- 
tions at x = +00 leading to degenerate kinks is given by 
H_/K. 

In addition to the degeneracy due to different bound- 
ary conditions, any kink solution will have an "internal" 
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symmetry group, denoted by /. This group will contain 
all those transformations that leave unchanged the whole 
kink solution (including the boundary conditions). So we 
have I C K and the space of degenerate kinks is H- / 1. 

The kink classification problem can be described in 
more detail as follows. Suppose we fix the boundary 
condition at x = — oo to be then the only con- 
straint on the boundary condition at x — +00 is that 
it should be in the distinct topological sector. There 
is a full vacuum manifold (mod Z2) worth of choices for 
$(x = +00) = $ + . For certain choices of <E>+ we can solve 
the equations of motion and obtain a set (described by 
the space K/I) of kink solutions that extremize the en- 
ergy. Let the value of this energy be U[$ + ; <&_] where we 
have explicitly indicated that different choices of bound- 
ary conditions can lead to kink solutions of different en- 
ergy. We are interested in the space of minima of the 
"potential" [/[<!>+; with respect to The global 
minima of this potential will describe the lowest energy 
kink solutions in the model and may be termed the "kink 
vacuum manifold" . Other local minima will describe kink 
solutions that are separated from the lightest kink by an 
energy barrier. In this way it might be possible to obtain 
"generations" of stable kinks having the same topological 
charge but differing in their energies. 

The existence of the kink solution found in the previ- 
ous sections does not tell anything about whether it is a 
minimum of {/[$+; $-]. To examine if the solution is a 
local minimum we can do a perturbative analysis where 
the perturbations are required to vanish at x — —00 - 
so as to hold $_ fixed - but are not required to van- 
ish at x = +00 - since we want to find changes in the 
energy when the boundary conditions are varied. Going 
back to Sec. [y] we find that all but one of the non-trivial 
Schrodinger potentials are positive (for N > 3) at both 
spatial infinities. This means that a perturbation that 
does not vanish at both spatial infinities, gives a diver- 
gent contribution to the energy. Hence the kink solution 
is stable under these perturbations of the boundary con- 
ditions. The only exceptional case is the off-diagonal 
perturbation of the type in eq. (|44|) . The potential (eq. 
(ff8|)) goes to a positive value at x — — 00 but vanishes 
from below at x = +00. In fact, as described below eq. 
( ff8| ) , there is a zero mode for this perturbation that van- 
ishes at x = —00 but goes to a non- vanishing constant at 
x = +00. So this mode is a "dangerous" one and needs 
to be examined further. 

A closer inspection of this mode shows that it corre- 
sponds to gauge rotations of the field $ which leave <&_ 
invariant but rotate 4> + . In other words, the zero mode 
rotates the kink within its own class described by the 
space H/I. Hence, the perturbation under consideration 
is not an instability but a gauge rotation. Therefore the 
class of kinks that has been found describes a set of local 
minima of $_]. 

There is a subtlety in the discussion above which we 



have glossed over. To obtain the Schrodinger equation in 
eq. ( |39| ) we have to perform an integration by parts and 
assume that the boundary contributions vanish. How- 
ever, here we are considering perturbations that do not 
vanish at infinity. This fact means that there is an extra 
contribution to SE^ in eq. (p9) given by 



SEoo = ~Tr 



dx 



- +00 



(55) 



The contribution at x = —00 vanishes because we are 
choosing ^(—00) = but the contribution at x — +00 
does not obviously vanish and depends on the derivative 
of the perturbation at infinity. However, since we are 
only interested in finite energy field configurations and 
the energy density contains the term Tr(9 3 ;5') 2 , we must 
require that the derivative of the perturbations vanish at 
spatial infinity. So it is rigorous to take SE^ = even if 
the perturbations do not vanish at infinity. 

A more extensive discussion of the kink classification 
problem is left for future work. 



VIII. SUMMARY 

We have considered SU (N) x Z 2 models. On the ba- 
sis of topology, the model will contain topological kink 
solutions. A general technique for constructing the solu- 
tions is not known. Here, by using some guesses and by 
choosing a special relation between parameters (eq. (0)), 
we have analytically constructed a class of kink solutions 
(see eq. (p8|)). The energy of the solutions is 



E = 



2V2 fN-l\ m 3 



3 A 



(56) 



The limiting value for large N is equal to the energy of 
a Z2 kink with mass parameter m and coupling constant 
A. We have explicitly checked that these kink solutions 
are perturbatively stable. It is not known if the solutions 
are globally stable and this remains an interesting open 
problem. 

We have then described the space of kinks as par- 
titioning into distinct classes. All members have the 
same topology, yet elements of different classes are not 
expected to have the same energy. The solutions con- 
structed above describe only one of the (unknown num- 
ber of) classes of kinks and might lie on the "kink vacuum 
manifold" - the manifold consisting of the least energetic 
kinks in the model. 

We hope that the solutions found here can be used as 
a guide to the construction of other topological defect 
solutions in complicated field theories. 
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